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Abstract 

The short-time evolution of a growing interface is studied analytically and numerically 
for the Kadar-Parisi- Zhang (KPZ) universality class. The scaling behavior of response and 
correlation functions is reminiscent of the "initial slip" behavior found in purely dissipative 
critical relaxation (model A). Unlike model A the initial slip exponent for the KPZ equation 
can be expressed by the dynamical exponent z. In 2+1 dimensions z is estimated from the 
short-time evolution of the correlation function for ballistic deposition and for the RSOS model. 

1. Introduction 

Interface formation and growth are typical processes in nonequilibrium systems. Two important 
examples are fluid flow in porous media Q and deposition of atoms during molecular beam epitaxy 
(MBE) B. It is expected that at times much later than typical aggregation times and on 
macroscopic length scales these interfaces develop a characteristic scaling behavior, where the 
scaling exponents fall into certain dynamic universality classes [Q, ||. In certain cases, however, 
interfaces can also show turbulent, i.e., spatial multiscaling behavior Usually a d-dimensional 
interface is embedded in d + 1-dimensional space such that the interface position at time t can 
be described by a height function h(x,t), where x denotes the lateral position in a ci-dimensional 
reference plane given by the surface of a substrate. Complete information about the scaling 
behavior is contained in the dynamic structure factor, which is related to the time displaced 
height-height correlation function C(x — x', t, t') = (/i(x, t)h(x' , t')) — (/i(x, t)) (/i(x', t')), where a 
laterally translational invariant system is assumed. For t,t' — > oo and finite \t — 1'\ the correlation 
function displays the asymptotic scaling behavior 

C(x - x', t, t') = |x - x'\ 2a F c (\t - i'|/|x - x'H, (1.1) 

where a denotes the roughness exponent and z is the dynamic exponent [|l], g]. For a laterally 
translational invariant system the interfacial width w 2 (t) = (/i 2 (x, t)) — (/i(x, t)) 2 is only a function 
of t and displays the scaling behavior w(t) ~ t@ for late times, where (5 = a/z is the growth 



exponent. For MBE as an example the scaling behavior displayed in Eq.(Ll) gives access to the 
exponents a and z both experimentally by reflection high energy electron diffraction (RHEED) 
(see, e.g., chaper 16 of Ref. |l]]) and by direct imaging using a surface tunneling microscope @ and 
theoretically by continuum models [[l], and Monte-Carlo simulations Q 

Continuum descriptions of interfacial growth processes can be obtained from general symmetry 
principles and conservation laws obeyed by the growth process For a wide class of growth 
processes the resulting continuum model is given by the well-known Kadar-Parisi-Zhang (KPZ) 
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equation 0, which reads 

f /i(x, t) = vV 2 h{x, t) + f (V/i(x, t)) 2 + r/(x, t). (1.2) 
The noise r/(x,t) has a Gaussian distribution with (rj(x,t)) = and 

(t/(x, t)r/(x', 0> = 2D5(x - x')<5(t - i'). (1-3) 
The parameters -D, and A are assumed to be constants and averages (...) are taken over the noise 



distribution. In the long time limit Eq.(L2) has a global symmetry which is commonly denoted as 
Galileian invariance |l], [jj. An important consequence is that the exponents z and a of the KPZ 
equation fulfill the scaling relation a + z = 2. The exponents of the KPZ equation are exactly 
known only in d = 1, where z = 3/2 and a = 1/2 due to the existence of a dissipation fluctuation 
theorem [J. In d = 2 numerical investigations indicate z ~ 1.6 and a ~ 0.4 [jjj. For d > 2 the 
asymptotic scaling behavior is either governed by linear theory (A = 0, weak coupling regime) or 
by another set of exponents inaccessible by analytical methods (strong coupling regime) depending 
on the value of the effective coupling constant g = DX 2 /(4i/ 3 ) jl], |j. Furthermore, it is interesting 
to note that the nonlinearity in Eq.( |1.2| ) renders all other possible nonlinearities irrelevant in the 
renormalization group sense in the long-time limit. For intermediate times, however, the presence 
of other nonlinearities in the growth equation gives rise to various crossover phenomena [|l|, |i~0| , 

2. Analytic Theory 



In Fourier space the KPZ equation Eqs.(L2) and (|l.3j ) is equivalent to the dynamic functional 



J = Jq + Ji H HI' 13 wn i cn consists of the Gaussian part 

Jo[h, h]=J J™ dt j£>fc(q, t)fc(-q, t) - fc(q, t) ( J^Hl, + ^^(-q, *)) } (2-1) 



and the interaction part 

A f d d qi f d d q 2 f°° 



Jl[h,h]=--J J ^-^^° ^q 1 -q 2 /i(-qi-q 2 ,t)/ i (qi,t)/i(q 2 ,t), (2.2) 



where h(q,t) is the Fourier transform of the response field |L2|. The initial condition /i(q, 0) = 0, 
which is implicitly assumed in Eqs. (|2.l[ ) and ( |2.2[ ), breaks the temporal translational invariance of 



the KPZ dynamics. The analytic treatment of Eqs^^l) and ( U-A ) is based on the identities 



-§- t h(q_,t = 0) = 2Dh(ci,t = 0) and G(q = 0, t, if < t) = 1, (2.3) 

where the response function G is given by the formal average {h{— q, t)h(q, t')) with respect to the 
dynamic functional J[h, h] p5| . For details of the field-theoretic treatment of Eqs.( |2.1| ) and ( p.2| ) 
we refer to Refs.|TT|, |l5| and only quote the main results for later reference. 

The correlation function C and the response function G are found to obey the scaling relations 
C(q,M' « t) = {t'/t) e \q\- d - 2a f c {\q\H) and G(q,M' « t) = {f /tf\q\- d f G {\q\H), (2.4) 
respectively, where in contrast to model A [p | the short-time exponents 9 and 6* are given by 

= {d + 4)/z - 2 = (d + 2a)/z and ^ = (2.5) 
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at the nontrivial fixed point (A 7^ 0) of Eq.(1.2). In d = 1 the exact value 9 = 4/3 can be obtained 
from the exact value z = 3/2 and is confirmed by a Monte-Carlo simulation of ballistic deposition 
[15]. From numerical estimates for z in d = 2 one obtains 8 ~ 1.7. The exponent relation given by 
Eq. Q2.5| ) simply means that the short-time and the long-time scaling behavior of the correlation 
function are identical, i.e., the short-time scaling behavior can be obtained by extrapolating the 
t'-dependence of C(q, t,t') from t' ~ t to t' = 0. In fact, the scaling relation given by Eq.(p~5|) 
can be derived independently by analyzing the fluctuation spectrum of the interface displacement 
velocity averaged over a macroscopic portion of the interfacial area [1£[. It should be noted, 
however, that the perturbative analysis of Ref . |i~5|l only consitutes a rigorous proof of Eq.( |2.5[) for 
d = 1. For d > 2 one encounters the strong coupling regime of Eq. (|1.2| ) which can no longer be 
treated analytically. 

Finally, we remark that an alternative scaling form for C can be obtained from the definition 
of the growth exponent f3 which leads to 6 = d/z + 2/5. The scaling behavior displayed in Eq.(|2.4|) 
can then be written in the simplified form C(q, t,t' <C t) = t'° gc(\({\ z t), where gc(y) = y~ e fc{y)- 



3. Monte-Carlo Results 

The scaling behavior of C(q, t, t' <C t) according to Eq.( |2.4| ) can be tested numerically by a Monte- 
Carlo simulation of simple deposition models on lattices The continuum description underlying 
Eqs.( |1.2[ ) and ( |1.3|) is replaced by a discretized description according to 

h(x,t) = h(x = aj,t = n/(FL d )) =ahj(n), (3.1) 

where the lattice constant a is assumed to be the same both in the plane of the substrate and 
perpendicular to it and j = (ji, . . . , jd)- The lattice has L d sites, F is the incoming particle flux, 
and n is the number of deposited particles. Furthermore, the incoming particle flux F has been 
normalized to unity, so that t in Eq.( |3.ip is dimensionless and given by the number of deposited 
layers. Finally, hj(n) defined by Eq. (|3.1| ) is also dimensionless and denotes the number of particles 
deposited at lattice site j after n particles have been deposited on the lattice. Ballistic deposition 
on a two-dimensional substrate is defined by the deterministic growth rule 

hj,k(n + 1) = max(/ij_i jjfc (n), h j<k ^ 1 (n), h jtk (n) + 1, h j+ltk (n), h jtk +i(n)), (3.2) 



(see, e.g., Ref. Q), where the site (j, k) in Eq. (|3.2[) has been selected randomly from the Lx L sites 
of the lattice and periodic boundary conditions are applied. 

In order to measure the short time exponent given by Eq.(|2.5|) it is sufficient to probe the inte- 



grated time displaced correlation function, i.e., one probes C(q = 0,t,t') as described in Ref.[15|. 
Like a real deposition process the simulation is characterized by an a priori unknown microscopic 
aggregation time t a . The scaling behavior of C according to Eq.( [2.4[ ) can only be observed for 
t' 3> t a . On the other hand t' <C t is required for Eq. fl2.4|) to hold, so that short-time scaling 
is restricted to the time window t a <^ t' t. Furthermore, the lattice size L must be chosen 
sufficiently large in order to avoid the onset of finite-size crossover effects if t' l l z ~ L when t' is 
still much smaller than t. For the simulation described here t = 1000 and L > 200 fulfill the above 
requirements. In order to cope with the very small signal to noise ratio in each measurement of 
C(0,t,i') for t' < t averages are taken over 2 x 10 4 realizations. These are distributed over 20 
individual runs at every point in time so that the jackknife method can be applied for the data 
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Fig. 1: Correlation function C(0,t,t') in d = 2 for ballistic deposition (top) and the RSOS 
model (bottom) as a function of t'/t for 0.005 < t'/t < 1 and L x L = 300 x 300 (solid line). 
Error bars are shown only at a few selected points in time and represent one standard deviation. 
The dashed lines display power laws with the measured short-time exponents 6 — 1.729 ± 0.042 
and 9 = 1.733 ± 0.041, respectively. The data follow this power law rather accurately in the 
interval 0.02 < t'/t < 0.2 (bold arrows). 



analysis. The result for ballistic deposition according to Eq. (|3,2| ) and for the RSOS model in 
d = 2 according to Ref . [[17|] is displayed in Figjl], where C(0,t,t') is shown as a function of t'/t for 
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t = 1000 and L = 300. According to FigJl] the interval 0.02 < t'/t < 0.2 is available to determine 
the short-time exponent 6. From Eq. fl2,5|) we obtain the estimate z = 1.608 ± 0.013 by averaging 
over the two values for 9 given in Figj|. Finally, we note that according to the scaling relation 
a + z = 2 one has a = 0.392 db 0.013 for the roughness exponent. These values are in agreement 
with other numerical data for z and a in d = 2 (see chapter 8 of Ref. Q for a collection of recent 
estimates) and they therefore provide some support for the general validity of Eq,( |2,5| ). 
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